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We give a brief review of the Kondo effect and exactly solve the nonequilibrium Kondo 
problem at the special point in the parameter space of the model with spin-dependent chem- 
ical potentials. Using the obtained solution, we compute several experimentally observable 
quantities: charge current, spin current and magnetic properties. Applications to the phys- 
ically important cases such as the case where the impurity spin is placed in between two 
ferromagnetically polarized leads are discussed. 
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1. Introduction 

In most metals, the resistivity monotonically decreases with the decrease of temperature 
because it is dominated by phonon scattering. On the other hand, the metals with magnetic 
impurities such as Fe have the resistance minimum at a certain temperature, known as Kondo 
temperature Tk, and under Tk the resistivity approaches unitary limit as the temperature 
approaches absolute zero. The experimental discovery of this effect dates back to the early 30s. 
The reason for the resistance minimum was not known at that time. In the 60s, however, a very 
significant advance of the theory for the effect was developed by Kondo. 1 Kondo's calculation 
of the resistivity, which was to explain this minimum, was based on the s-d model where the 
local impurity spin S is coupled via an exchange interaction J with the conduction electrons 
of the host metal. Using third order perurbation theory in coupling J, he showed that this 
s-d interaction leads to singular scattering of the conduction electrons near the Fermi surface 
and a InT contribution to the resistivity. The InT term increases at low temperatures for an 
antiferromagnetic coupling and when this term is included with the phonon contribution to 
the resistivity, i.e., T 5 term, it is sufficient to explain the observed resistance minimum. Since 
his perturbation calculations could not be valid at low temperatures, a more comprehensive 
theory was needed to explain the low temperature behavior. This task attracted a lot of theo- 
retical interests in the late 60s and early 70s. More extensive perturbative calculations which 
sum the most divergent terms were performed by Abrikosov. 2 The scaling idea introduced by 
Anderson provided a new theoretical framework for the Kondo problem. 3 This idea was taken 
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over by Wilson's renormalization group method discussed below. Another perturbative ap- 
proach was developed in the study of the Anderson model, which is more microscopic than s-d 
model, where perturbative calculations for Coulomb repulsion U have no divergence. 4 Nonper- 
turbative approaches to the Kondo problem will be overviewed below. From the present point 
of view, the Kondo effect is understood as a precursory phenomenon of the formation of the 
bound state between the conduction electrons and the impurity spin, known as Kondo singlet. 
The theoretical structure underlying the Kondo problem is similar to the one in Quantum 
Chromodynamics (QCD) and the Kondo effect corresponds to the confinement of quarks and 
gluons. 

Recently, the Kondo effect is revived in the context of mesoscopic physics such as semi- 
conductor quantum dot devices where the quantum dot spin and the electrons in electrodes 
corresnpond to the impurity spin and the conduction electrons in the original Kondo problem, 
respectively. 5 The crucial difference between the usual Kondo effect and that in mesoscopic 
systems is the nonequilibrium effect, i.e., the voltage vias between electrodes. It is very diffi- 
cult to approach theoretically to the nonequilibrium situation of the Kondo problem. One of 
the standard method to treat nonequilibrium problems perturbatively is the Keldysh Green's 
function method (discussed in ref. 6 in connection with the Kondo effect). Nonperturbative 
approaches will be discussed in detail below. 

Before discussing the nonequilibrium Kondo problem, let us briefly overview the nonperturba- 
tive approaches to the equilibrium Kondo problem (more detailed information can be found in 
ref. 7 and 8). The first important contribution developed by Wilson in the 70s is the Numerical 
Renormalization Group (NRG) method 9 which provides a way to numerically perform scale 
transformations proposed by Anderson. By using NRG, Wilson has obtained definitive results 
for the ground state and low temperature behavior for the spin half s-d model. The second 
development in the late 70s and early 80s is the discovery of exact solutions to the Kondo prob- 
lem. By using Bethe Ansatz (BA) method which is a powerful tool to solve one-dimensional 
integrable models such as the Heisenberg model and the Hubbard model, Andrei 10 and Wieg- 
mann 11 solved the spin half s-d model independently of each other. The critical idea, already 
used in Wilson's NRG, in BA approach is the mapping to the equivalent one-dimensional 
model. After these works, Wiegmann showed that BA approach could be extended to give ex- 
act results for the Anderson model. 12 The most recent nonperturbative approaches by Affleck 
and Ludwig in the 90s are based on Conformal Field Theory (CFT). 13 They have formulated 
the Kondo problem as a boundary quantum critical phenomenon and have applied CFT to 
it. First, they found the nontrivial strong-coupling fixed point where the absorption of the 
impurity spin is elegantly described by the deformation of Kac-Moody algebra, infinite di- 
mensional Lie algebra, for the spin part. Second, they assumed that the absorption process at 
the point away from the fixed point is governed by Kac-Moody fusion rules, the analogue of 
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addition of angular momentum in CFT, and extracted the anomalous non-Fermi liquid(NFL) 
behavior of the multichannel Kondo problem introduced by A. Blandin and P. Nozieres. After 
the works of Affleck and Ludwig, Emery and Kivelson have mapped the two-channel Kondo 
problem into the resonant-level model, which is equivalent to noninteracting fermions for a 
particular value of the z component of the exchange similar to the Toulouse limit, 14 by using 
the standard Abelian bosonization technique. 15 Since the resultant resonant-level model is a 
quadratic form in electrons and impurity spin operators at this special Fermi Liquid (FL) 
fixed point, the absorption of the impurity spin is described by the unitary transformation 
which diagonalizes the above-mentioned quadratic form. 

We shall now leave the equilibrium case and turn to the nonperturbative approaches to the 
nonequilibrium Kondo problem. The first contribution is given by Schiller and Hershfield. 16,17 
They combined the y-operator method 18 developed by Hershfield with the Emery-Kivelson 
solution discussed above and then obtained the exact solution for the nonequilibrium Kondo 
problem at a special point in the parameter space. After this work, Delft et al have de- 
veloped the CFT approach combined with Hershfield's y-operator formulation. 19 Recently, 
Mehta and Andrei have applied the Nonequilibrium Bethe Ansatz (NEBA) to the computation 
for steady state properties of quantum impurity problems such as interacting resonant level 
model(IRLM). 20 Here we should note that the comprehensive nonperturbative understanding 
of the nonequilibrium Kondo problem is still an open problem and requires more extensive 
studies. 

In this article, we extend the Schiller-Hershfield solution to the nonequilibrium Kondo 
problem to include the case where the chemical potential of each lead depends on the spin 
of conduction electrons. Actually, these situations are often discussed in the area such as 
spintronics. It is different from Ref.16 and 17 that we employ a description of nonequilibrium 
steady states (NESS), 21 ' 22 developed in the context of mathematically rigorous theory such 
as C*-algebra, instead of using Hershfield's y-operator. The obtained results given below are 
generalization of those in ref. 16 and 17 and reproduce their results as limiting cases. 

2. Model and mapping 

The physical system under consideration is shown schematically in Fig.(l). It consists of 
left(L) and right(R) leads of noninteracting spin-1/2 electrons, which interact via an exchange 
coupling with a spin-1/2 impurity spin placed in between the two leads. Following Ref., 13 
the conduction-electron channels that couple to the impurity are reduced to one-dimensional 
fields ip a a(x), where a = L/R and a =|, [ are the lead and spin indices respectively. In what 
follows, we assume that the conduction-electron dispersion around the Fermi level is linear, 
i.e., efc = hvpk, where and k are measured relative to the Fermi level and Fermi wave 
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Fig. 1. Schematic description of our system. A tunnel junction consists of two leads of noninteracting 
spin-1/2 electrons and a spin-1/2 impurity spin placed in between the two leads. The energy levels 
denote the spin-dependent chemical potentials of the conduction electrons. 

number respectively. The most general form of the s-d Hamiltonian H is given by 

lpia( X )-Q-*Paa(x)dx 
a=±j,rt a= | ,4, D X 

+ E E jf^- mgi HT\ (1) 

a,/3=L,R X=x,y,z 

where f denotes the impurity spin and s^g is defined as 



S X a(3 = T,^A0)^f^A0), (2) 

a, a' 

where a x , (A = x,y,z) are a Pauli matrices. And \xb and gi are the Bohr magneton and the 
impurity Lande g factor, respectively. In addition to the above Hamiltonian H, we consider 
the following nonequilibrium conditions: 

Y = Y e + Y s , 

^ roc 

Y = ~T~ E / WhJta ~ lj}^Ra\dx, 

Y = —Y^^La^La-V^RatpR^dx, (3) 
cr 

where V e and V s correspond to a usual electronic voltage and a spin voltage respectively. 
In Ref. 16 and 17, Schiller and Hershfield have considered the case where only V e exists. On 
the other hand, there are some situations, such as ferromagnetic leads cases, where the spin 
voltage V, exists and plays a key role in the spin transport phenomena such as spincurrent. 
Therefore, the nonequilibrium Kondo problem with the above more general nonequilibrium 
conditions needs to be newly formulated. 

We shall discuss how to diagonalize the above Hamiltonian. First, we restrict the region of the 
jf parameter space to the case where J£ a = J£ a = Jf, j£ R J^ R = J^ R = J RL = 0, 
and J^ L = J RR = J z takes a particular value (discussed below) . We note that setting J^ R = 
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produces no significant difference on the low-energy physics since, in the renormalization 
group sense, J% R 7^ is generated on scaling from weak coupling. 23 Next, we bosonize eq.(l) 
and eq.(3) respectively by following the Emery-Kivelson solution of the two-channel Kondo 
problem. 15 We introduce four different boson fields to account for four different fermion fields 

as 

V 2ira 

= -^(-l)%+% e -^W, (5) 

V lira, 

^ t (x)~(-l)%e- !$ «tW, (6) 

\/2TTa 

^(x)~e- !$ «»W, (7) 
V 2na 

where a is a high-energy cut-off (lattice spacing), and number operators for fermion 7Vj )(T (i = 
L, R, a =|, I) are defined as 

J dx : Vl(x)Vw(x) := — J dx—Q—, (8) 

and satisfy the canonical commutation relations [$>i a (x), Nj a ] = iSij5 acr i. The boson fields $j jCr 
also satisfy standarad commutation relations 

[®ia{x), ®ja'(y)} = -iTTSijS aa/ sgn(x - y), (9) 

and hence we can make sure of the following correspondence between the fermion fields and 
boson fields: 

27 r : liWMz) : = (10) 

ivf J_Jx : lU*)^lM*) : = -J_Jx: (-^— ) : (11) 

where : AB... : denotes the normal ordering of operator products AB.... Furthermore, we 
introduce four new boson fields constructed by the above four boson fields in eq.(7) as 

$ c = + ®li + §m + ®Rl)i 

®s = - ®li + $m - ®w)> 

= \(®L] + ®Ll ~ <$>R] - &m), 
®sf = -($L T -$L|-$ilT + $ill), ( 12 ) 

where subscripts c, s, f, sf correspond to collective charge, spin, flavor (left minus right), and 
spin-flavor modes respectively. Now, the Hamiltonian and nonequilibrium conditions can be 
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H 



flVF \ -> f 

~4nF ^ J 



dx 



v=c,sj,s}' 



dx 



+ 



J+ 
ira 
J- 
ira 1 

Tra 



[-T x sm(xs) + t s/ cos(x s )]cos(x s/ ) 
[t x cos(x s ) + T y sm(xs)]sm(xsf) 
t x cos(xs) + T y sm(x s )}sm(xf) 



(13) 



and 



Y e = 
Y s = 



eV e r°° d& f {x) 
2vr dx 

eV s f°° , d$ a f(x) 



2tt /_ 



dx 



dx 



(14) 



where are the even and odd combinations of J± L and J RR , i.e., J ± = (</f L + J RR )/2 and 
XsiXf an d Xsf ar e defined as 

Xs = $ s (0)-|(7V c -AT s ) + |, 

7T 



Xf = ^f(0)-^(2N C -N f -N sf ), 



7T 



Xsf = ^sf(0)--(N f -N sf ). 



(15) 



The next step is to employ the canonical transformation H' = UHU^ and Y' = UYqU^ with 
U = exp(ix s r z ). Since Xs commutes with both Xf an d Xsf, the transformed Hamiltonian H' 
is written as 



v=c,s,f,sf 
jLR 

Tra 



J+ 



J- 



+ — T v cos(x s f) T x sin(x s /) 

\ ox J Txa 1 ixa 



d$ s (0) 



dx 



where we have used the commutation relation 

3* 



-27ri^5(a;)e iXs 



(16) 



(17) 



and have neglected divergent constant (the validity has been discussed in ref. 24). The nonequi- 
librium conditions Y e and Yq are unaffected by the above canonical transformation, i.e., 
Y' = Yq. The last step is refermionization of the above bosonized operators. To this end, 
we first express the spin f in terms of fermion operator as d = ir + . The four boson fields 
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<& u {y = c, s, f, sf) are refermionized according to 



' -i(*/(x)-f (2N C -N f -N sf )) 



4> f (x) = _ e -n-/w-, 

V lira 



^ sf{x) = « e -i(*./(*)-5 (18) 
V27ra 

etc. Here the phase iird)d, similar to the Jordan- Wigner transformation, guarantees the anti- 
commutation relations between d and ip u (x). 

Once these steps are completed, we can obtain the refermionized forms for H' and Y': 

J+ 



/oo a j+ 

dx ^l {x) Mx) + — f (0) + ^/(0)](dt - d) 

+ -±=[$(o) - V/(o)](dt + d) + -=[^t (o) - ^(o)](dt + d) 

2\/2ira 1 2\l2ixa 1 

+ [^tf - (J 2 - 2tt^ f ) : Vj(O)VaO) :] (Vd - ^ , (19) 



and 



/oo 
-oo 

/oo 
dx*l>\ f {x)ij> sf {x). (20) 
-oo 

Here we can immediately notice that (19) would be a quadratic form in fermion operators 
if we set J z = 2ttHvf- Since the nonequilibrium condition Y' is originally a quadratic form, 
both H' and Y' reduce to quadratic form at this special point. Hence the strongly interacting 
nonequilibrium problem maps to a noninteracting one. For convenience, we introduce new 
Majorana fermions as 

d + d) t rft - d 

a = ^r- "=ivt- (21) 

Here we note that Majorana fermions satify d 2 = b = 1/2. After Fourier transformation, H' 
and Y for J z = 2-Khvp can be rewritten as 

jLR j- 

where L is the size of the system, = ftvpfc and we have neglected the terms involve the ip c 
and tp s which are decoupled from d and b. 
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3. Solution and observables 

In this section, we shall show how to diagonalize eqs.(22) and (23) simultaneously. After 
obtaining the exact solution, we will calculate experimentally observable quantities such as 
charge current, spin current and magnetic properties. 

Since the precise procedure for diagonalization of H' is given in ref. 17, we shall only outline 
the main steps briefly. The first step is to construct scattering states c uk which diagonalize 
H' and are defined by the operator equation 

ici, k ,H'} = -e k c\ k + itjtyt fc _ c l k ) [y = /, 8f)t (24) 

where the positive infinitesimal rj is introduced to guarantee appropriate boundary conditions. 
Detailed derivation of the scattering-states operators are given in the APPENDIX B in ref. 
17. Using c U) k, our Hamiltonian H' can be written in a diagonal form as 

H ' = E ( 25 ) 

Next we need to expand original operators, ip^l, a and b, in terms of scattering-states opera- 
tors c^\. TABLE II in ref. 17 summarizes the results of expansion. At this point, we employ a 
description of nonequilibrium steady states (NESS) instead of the ^-operator method devel- 
oped by Hershfield. In NESS approach, once we can obtain steady state density matrix p + , 
we impose the following conditions on the two-point function: 

i c l,k c ^k') = F^e^S^Skw (n, v = f, sf), (26) 

where (...) = Tr(...p + ) stands for the NESS average and F M (e) are defined by the ordinary 
Fermi distribution function /(e) as Ff(e) = f(e — eV e ) and F s f(e) = f(e — eV s ) respectively. 
We should note here that there is a somewhat subtle point to use C*-algebraic approach, since 
we involved non-local operators such as iVj in our bosonization procedure. 

Now, we are ready to calculate the observable quantities. The physical observables of 
interest such as charge current I c , spin current I s and impurity magnetization M z are rewritten 
after the canonical transformation U : 

jLR 

= B*^ + '^ )& ~ - *■>**• ( 28 > 

M z = iuBQiab, (29) 

(the derivation can be found in ref. 17). From now on, we focus on the average of transformed 
operators, since (...) = lim t _ >+00 Tr(... e - im p e iHt ) = lim t ^ +0O Tr(C/...L7t . Ue~ iHt U^ ■ Up W ■ 
Ue iHt tf) = \im t ^ +oo Tr(U...Uh~ iH,t p' e iH,t ) = (U...U*), where p is the density matrix for 
the initial state and we have used the cyclic property of the trace. Using the condition for 
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To 

1 2 


( J-) 2 lA-KnfiVT? 

\tj 1 i t: l \ Li/I LU J-t 


r L 




r« 


{jf R ) 2 /ATtahv F 


J- rrt 


J± R J~ /4irahv F 




jf R J+ /4nahv F 



Table I. Symbols and their definitions. 



the two-point correlation function (26), we can obtain the average of I' C ,I' S and M z due to 
the quadratic nature of them. (Generally speaking, we can also compute higher-order terms 
by using Wick's theorem for c v ^). 

Charge Current and Spin Current 

First, let us discuss the charge current and spin current in our system. Since the charge 
current obtained here is exactly equal to the result in, 17 we only write down the result: 

I C (V) = (I' C (V)) = J ^ deA a (e)[f(e - eV e ) - f(e + eV e )}, (30) 

where 

Me) = - Im {( e + ^ (e +trl)- ( ^) 2 } ' (31) 

and T a , (a = l,a,b) are tabularized below. Here, we note that I C {V) depends only on V e and 
is independent of the spin voltage V s . 



We will now leave the discussion of charge current and turn to that of spin current. 
We obtain the following result which includes the one in ref. 17 as a specific case: 

Is(V) = (I' s ) 



1 



2irh J_ 



de 



f(e - eV e ) + f(e + eV e ) f(e - eV s ) + f(e + eV s ) 



(ji B giH)(T L - T R )Tie 



\(e + ir a )(e + ir b )-( PBgi Hy\ 2 

+ ^ de[f(e - eV s ) - f(e + eV s )\ \T 2 A a (e) + T b A b (e)] 
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Fig. 2. The spin differential conductance G S (V S ,T) as a function of V (left) for zero magnetic field 
and different temperatures and (right) at T — and nsgiH — 2T a} for different ratios of T a to 

rv 



d e[f(e-eV s )-f(e + eV s )} 



2(e 2 + T a r b + (p b9l H 



27rhJ_ 00 uy ' JK n |(e + ir a )( e + ir fe )-(^i/) 2 | 2 ' 

(32) 

where A b is defined as 

Me) = " Im { (e + i r a )(e^Z a -(»B 9 .Hy } ■ (33) 

We can immediately notice that eq.(32) reproduces eq.(6.7) in ref. 17 if we set V s = 0. We 
also notice that we have finite spin current I s even though the case where J^ L = J RR , i.e., 
symmetric case due to the existence of IV 

Let us now discuss the spin current for specific cases, (i)symmetric, zero magnetic field, at 
finite temperature and (ii)symmetric, finite magnetic field, at T = in full detail. If we set 
jf L = J RR , Tl = Fr and 1^ = 0, then we can write down I s as 

Is(V s ) = ^ j deA b (e)[f(e - eV s s) - /(e + eV s )]. (34) 

Contrast to the charge current case (30), I s involves A b , i.e., the spectral function for the 
Majorana fermion b . We can obtain a closed form expression for (34) using the digamma 
function ^(z) 2 ^ which is defined by the gamma function T(z) as 

1>(z) = inr(z). (35) 

For example, case (i), the spin current is rewritten as 

^ImWU^n, (36) 



nh { \2 2-nk B T 
where ks is Boltzmann constant. For finite magnetic field case, we can also explicitly rewrite 
I s using the digamma function, however it takes a more complicated form. 
Next, we will show the differential conductance for the spin current G S (V S ,T) = dI s /dV s as 
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a function of bias, for (i) and (ii). For (i), the result is 

2e r ±_ B J,m(L + lL +ieV ' 



G ' (v "' T > = T2^ Re {#" (2 + "W" ) ) • < 37 > 

where ip^-\z) = di(j(z)/dz denotes the trigamma function. In the left figure of Fig. (2), we show 
G S (V S , T) for different temperatures. In the case (ii), we can easily carry out the integral, since 
the Fermi distribution function becomes the step function at zero temperature. Hence, the 
result is given by 

2e / eV a + iV a 

h brn \(eV s + iT a )(eV s + iT b )-^ Bgi Hf 

In the right figure of Fig. (2), we show G s (V s ,0) for a moderately large magnetic field 

fiB9iH = 2Tb at different ratios of r a to T b . For r a = r&, magnetic field splits the zero-field 

resonance into two Lorentzians centered about ±fis9iH. 



G S (V°, 0) = --r 6 Im ( ^,^^7^ t—— 2 ) . (38) 



Impurity Magnetization 

Next, we shall briefly discuss the impurity magnetization which provides direct infor- 
mation about the onset of Kondo screening. The average of the impurity magnetization 
(M z ) = {iuBQidb) can be computed as a function of the magnetic field H, the voltage V e 
and the spin voltage V s . The resultant expression is 



M(H, V e ,V s ) 

(r„ + r b )e 



a + Tb \ vr J_ 00 



r a + r 6 V vr |( e + ir a )(e + ir fe )-(^ J ff) 2 | 2 



f(e-eV" 



Tg + T b - £1 / LiB9iH (Tq + T b )e \ _ 

T a + T h V * e |( e + ir a )(e + ir fe )-(^^)2|2;^ e 6 > 

1 r L -r fi e 2 + r a r fe + ) 2 

+ I. de |( e + ,r a ) (e + ,r 6 )-(^) 2 | 2[/(e " eV > " /(e + ey W- 

(39) 

Setting y s = 0, eq.(39) reproduces eq.(8.2) in ref. 17. We can also extract more explicit 
expression from eq.(39) using the digamma function ip(z). 

4. Conclusion 

In conclusion, we have exactly solved the nonequilibrium Kondo problem at the special 
point in the parameter space of the model with the spin-dependent chemical potentials. Using 
this solution, we have computed several experimentally observable quantities: charge current, 
spin current and magnetic properties. We have found that the spin current does not vanish 
even in the case J± L = J± R , i.e., perpendicular part of the exchange interactions between 
two leads and spin are symmetric. Moreover, we have explicitly calculated the observables for 
the specific case in which the impurity spin is placed between two ferromagnetically polarized 
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leads. This case, a common situation in spin transport phenomena, has never exactly been 
discussed out of equilibrium and under the Kondo temperature before. In the future, it will 
be interesting to study the nonequilibrium Kondo effect in a quantum dot coupled to two 
noncollinear ferromagnetic leads. 26 
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